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Abstract 
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<^ Using the nonrelativistic QCD factorization formalism, we calculate the 

°" 

electroproduction cross sections of r\ c in ep collisions, including the contribu- 

Q_i| tion from both the transverse photon and the longitudinal photon. For this 

Oh' 

<]} [ process the color-singlet contribution vanishes up to the next to leading order 



perturbative QCD calculations. The dominant contribution comes from the 



color-octet subprocess. The nonperturbative color-octet matrix element 

of 1 5g 8 ^ of 7] c is related to that of 3 s[ 8 ^ of J ftp by heavy quark spin symmetry, 

and the latter can be determined from the direct production of J/tp at large 

transverse momentum at the Fermilab Tevatron. The measurement of this 

process at DESY HERA can be viewed as another independent test for the 

color-octet production mechanism. 
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I. INTRODUCTION 



Studies of heavy quarkonium production in high energy collisions may provide important 
information on both perturbative and nonperturbative QCD. Recent progress in this area 
was stimulated by the experimental results of the Collider Detector at Fermilab (CDF) 
at the Fermilab Tevatron. During the Tevatron run I, the CDF data [|TJ for the prompt 
production rates of J/ ip and ip' at large transverse momentum were observed to be orders of 
magnitude larger than the lowest order perturbative calculation based on the color-singlet 
model |J. At the same time, a new framework in treating quarkonium production and 
decays has been advocated by Bodwin, Braaten, and Lepage in the context of nonrelativis- 
tic quantum chromodynamics (NRQCD) |y. In this approach, the production process is 
factorized into short and long distance parts, while the latter is associated with the nonper- 
turbative matrix elements of four-fermion operators. So, for heavy quarkonium production, 
the quark-antiquark pair does not need to be in the color-singlet state in the short distance 
production stage, which occurs at the scale of 1/rriQ (rriQ being the heavy quark mass). At 
this stage, it allows a color configuration other than the singlet for the heavy quark pair, that 
is the color-octet. The later situation for heavy quarkonium production is called the color- 
octet mechanism. In this production mechanism, a heavy quark-antiquark pair is produced 
at short distances in a color-octet state, and then hadronizes into a final state quarkonium 
(the physical state) nonperturbatively. The color-octet terms in the gluon fragmentation to 
J/i()(ip') have been considered to explain the J/i()(ijj') surplus problems discovered by CDF 
fH|. Taking the nonperturbative matrix elements (0^^( 3 Si)) and (Of ( 3 Si)) as input pa- 
rameters, the CDF surplus problem for J/ip and ip' can be explained as the contributions 
of color-octet terms due to gluon fragmentation. In the past few years, applications of the 
NRQCD factorization formalism to J/ip(ip') production at various experimental facilities 
have been studied. 

Even though the color-octet mechanism has achieved some successes in describing the 
production and decay of heavy quarkonia, more tests of this mechanism are still needed. 
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Recently, the photoproduction data from the DESY ep collider HERA |J posed a ques- 
tion about the color-octet predictions for the inelastic photoproduction of J ftp More 
recently, possible solutions for this problem have been suggested in p|-|TTf. In this situa- 
tion it is certainly helpful to find other processes to test the color-octet mechanism in heavy 
quarkonium production. Under this context, we have studied r\ c photoproduction in NRQCD 
We find that this process has several unique properties for the test of the color-octet 



production mechanism: (1) this process is purely a color-octet process up to the next to 
leading order in perturbative QCD; (2) this process is dominated by the Sq channel, for 
which the color-octet matrix element can be related to the J/ip hadropro duct ion at large px 
by the heavy quark spin symmetry; (3) unlike the inelastic J/ip photoproduction which is 
affected by large diffraction background, the inelastic photoproduction of r] c has much lower 
diffraction background. So, the measurement of this process may clarify the existing conflict 
between the color-octet prediction and the experimental result on the J/ip photoprodution. 
In this paper, we will extend our previous study of rj c production from the photoproduction 
region {Q 2 rs 0) to the deep inelastic scattering (DIS) region (Q 2 > 0) at the energies rele- 
vant to the HERA collider. The electroproduction of r] c has the same features as the above 
mentioned for the test of the color-octet production mechanism. 

Moreover, since Q 2 can be large, electroproduction is a better process from which to 
test the color-octet mechanism and to extract the NRQCD long distance matrix elements 
than photoproduction. The latter process lacks any large scale other than the charm quark 
mass, and consequently, perturbative corrections to leading order calculations are large. In 
addition, nonperturbative effects, such as higher twist corrections to the parton model, are 
less effectively suppressed in photoproduction than in electroproduction at large Q 2 . 

The rest of the paper is organized as follows. In Sec. II, we give the cross section 
formula of rj c electroproduction in the color-octet model. In Sec. Ill, we give the numerical 
results about the production rates of rj c at the energy range relevant to HERA collider. The 
summary and discussions are given in Sec. IV. 
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II. rjc ELECTROPRODUCTION IN THE COLOR OCTET MODEL 



The power of the NRQCD formalism stems from the fact that factorization formulas 
for observables are expansions in the small parameter v , where v is the average relative 
velocity of the heavy quark and antiquark in quarkonium bound state. For charmonium 
v 2 ~ 0.3 , and for bottomonium v 2 ~ 0.1 . NRQCD velocity-scaling rules JTBJ allow us to 
estimate the relative size of various NRQCD matrix elements. This information, along with 
the dependence of the short-distance coefficients on a s and a, permits us to decide which 
terms must be retained in expressions for observables to reach a given level of accuracy. At 
low orders, factorization formulas involve only a few matrix elements, so several observables 
can be related by a small of parameters. 

In NRQCD the Fock state expansion for r) c is 

\rj c ) = O(l)\ccCS ,l)) + 0(v)\cc( 1 P 1 ,S)g) + 0(v 2 )\cc( 3 S 1: S or l)gg) 

+O(v 2 )\cc( 1 S ,8orl)gg) + ---. (1) 

For the production of r] c , the contributions to the NRQCD matrix elements from the last 
three terms of the above expansion are the same order of v 2 according to the NRQCD 
velocity scaling rules. They are all suppressed by v 4 compared to the contribution from 
the first term (the color-singlet contribution). However, the color-singlet contributions to i] c 
production vanish in the leading order and the next to leading order ^*g fusion processes 
(see the following). So, the production of r\ c is purely a color-octet process, even to the 
next to leading order of QCD calculations. Therefore, the r\ c production processes such 
as at the HERA, will provide an important test for the color-octet production mechanism. 
Furthermore, we will show by the following calculations, the dominant contribution to r] c 
production comes from the last term (color octet) of the Fock state expansion in Eq. (|1]). For 
this term, the associated color-octet production matrix element {Og c ( l So)) can be related 
to the matrix element (Cg ( 3 «Si)) by the heavy quark symmetry. And the latter color-octet 
matrix element is important for the color-octet explanation of the prompt J/ijj surplus 
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production at the Tevatron @,[|. So, the measurement of r\ c production processes is closely 
associated with the test of the color-octet gluon fragmentation mechanism proposed in 
The complete ep —>■ r] c X cross secion can be written as 



d<y ep ^ Vc x{s) 
dydQ 2 



T(Q 2 ,y)da^ (W 2 ). (2) 



Here T(Q 2 ,y) is the 



and 



°{l*P ^) = <r T (7*P + ^(7*P -> (4) 

where o"t and cr^ are the cross sections of virtual photon processes for the transverse polarized 
photon and longitudinal polarized photon respectively. 

According to the NRQCD factorization formalism ||, the production process 7* + g —>■ 
rj c + X can be expressed as the following form, 

da(^ + g ^r lc + X)=Y J F{l* + 9^n + X)(0%). (5) 

n 

Here, n denotes the cc pair configuration in the expansion terms of Eq. (HD (including 
angular momentum 2S+1 Lj and color index 1 or 8). -F(7* + g —>■ n + X) is the short distance 
coefficient for the subprocess 7* + g —>■ n + X. (O 1 ^) is the long distance non-perturbative 
matrix element which represents the probability of the cc pair in n configuration evolving 
into the physical state t] c . The short distance coefficient F can be calculated by using 
perturbative QCD in expansion of powers of a s . The long distance matrix elements are 
still not available from the first principles at present. However, the relative importance of 
the contributions from different terms in Eq. (|5|) can be estimated by using the NRQCD 
velocity scaling rules. 

From Eq. ([I]), we can see that the color-singlet contribution to the production of i] c is at 
the leading order in v 2 . The associated short distance coefficient is given by the subprocess 



1 *g^ccCSo,l)+9. (6) 

This process occurs at the next to leading order in a s for the j*g fusion processes. However, 
there is no contribution from this process, because it violates C (charge) parity conservation. 
(The C parities of the two gluon system (constrained in color- singlet) and the cc pair in 
( 1 5'o, 1) state are both +1, while the C parity of the photon is —1). 

The color-octet contributions to the r\ c production come from the leading order and the 
next to leading order 7* g fusion processes as shown in Fig.l. At the leading order in a s , the 
subprocess is 2 — ^ 1 (Fig. 1(a)), 

7 *<7-cc( 1 So,3). (7) 

For this process, we readily have 

■h 2 

atf + g - cc(% 8) - tfc) = ^-^S(s - 4m 2 - Q 2 )^ 1 ^)), (8) 

where s is the invariant mass squared of the partonic process. m c is the charm quark mass, 
and we have approximated the charmonium bound state mass of t] c by 2m c . 

At the next to leading order in a s , the subprocesses are 2 —>■ 2 (Fig. 1(b) and Fig. 1(c)). 
These processes contain the contributions from 

7* +g ->ccCSo,8)+g, 
7* + g -> ccCS^ 8) + g, 
1 *+g-*cc( l P 1 ,8)+g. 

The cross sections for these 2 — > 2 subprocesses can be expressed as the following form, 

+g^Vc + X) = ^FC s ^Lf)(Or( 2S+1 Lj)), (12) 
dt I671V 

where t — (z — l)s, and z is defined as z = p ■ k Vc /p ■ k* with p, k Vc , fc 7 « being the momenta 
of the proton, the outgoing r] c and the incident photon respectively. 

It is convenient to use the helicity amplitude method to calculate the cross section for- 
mulas for the virtual photon processes Eqs.(Q-(TT]). To the polarization vectors of the virtual 
photon, we choose 
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(9) 
(10) 

(11) 



1 1 Q 

&r= -t=(0,1,±i,0), e L = — pi + p 2) (13) 

V 2 Q Pi- P2 

for the transverse and longitudinal polarized photons respectively. For the incident and the 



outgoing gluons, following |L4j] , we choose their polarization vectors as 

(±) = N e [fa fa 4(1 T 7s)+ 4 fa fail ± 7s)], (14) 
{±) = N e [fa 4 fa(l t 7s) + fa 4Mi± 7s)]. (15) 

Where q = Pi + 2 ®. P2 P2 > Pi = ~Q 2 an< ^ P 1 * s ^ ne momenta for the incident photon. £>2, Ps 
and e2, e3 are the momenta and the polarization vectors for the incident gluon and outgoing 
gluon respectively. The normalization factor N e is 

= , 1 =, (16) 

V /2(Q 2 M 2 + su)t 

where M is the mass of rj c , the Mandelstam invariants s, t, u are defined as 

s = (pi+P2) 2 , i=(p 2 -p 3 ) 2 , u=(pi~P3) 2 , (17) 
and they satisfy the relation 

s + t + u = M 2 - Q 2 = Am 2 c ~ Q 2 - 
With these definitions for the polarization vectors of the photon and gluons Eqs. (|13| , [TJ], 



|15l) , the calculations of the helicity amplitudes are straightforward. The expressions for Fs 
in Eq.(|12|) are given in the Appendix. 



III. NUMERICAL RESULTS 

For the numerical evaluation, we choose m c = 1.5 GeV, and set the renormalization scale 
and the factorization scale both equal to /i 2 = (2m c ) 2 + Q 2 . For the parton distribution func- 
tion of the proton, we use the Gliick-Reya-Vogt (GRV) leading order (LO) parametrization 
0- 
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Because there is no color-singlet contributions to r\ c electroproduction up to the next to 
leading order perturbative calculations, we present in the following the numerical results 
only coming from the color-octet contributions. For the color-octet r\ c production, by us- 
ing the heavy quark spin symmetry, we estimate the associated color-octet matrix element 
(OffSo)) to be 

(OfCSo)} « (0?( 3 Si)) = 1.06 x 10- 2 GeV 3 . (18) 

The value of the color-octet matrix element (Og^ ( 3 Si)) follows the fitted value in [[HJ by 
comparing the theoretical prediction of direct J/ip production to the experimental data at 
the Tevatron. For the other two color-octet matrix elements, we estimate their values by 
using the naive NRQCD velocity scaling rules, 

(On'S,)) » (OTCSo)), (19) 

CTCpO) « (ofCs )). (20) 

In Fig. 2, We first display the contribtuion from 2 — > 1 partonic subprocess to the elec- 
troproduction of 77c, in which the intermediate color-octet state is 1 Sq 8 \ The 2 — ► 1 process 
contributes to the r] c production in the forward direction (z ~ 1). In e + p — > e + r] c + X, 
we study the inelastic rj c production in the kinematic range of 4 < Q 2 < 80GeV 2 and 
40 < W < 180Ge^ 2 for HERA experiments. 

Next we study the 2 — * 2 subprocesses contributions to the electroproduction of r) c . 
First, we consider the virtual photon-proton processes: ^*p — > T] C X. In Fig. 3, we plot 
the z distributions of the cross sections cr to t{l*p) = <^t(i*p) + °"-l(7*p) °f Vc production 
in the virtual photon-proton collisions at different values of Q 2 . Relevant to the HERA 
energy range, we choose the photon-proton cm. energy W = lOOGeV and select Q 2 = 
4, 10, AOGeV 2 to see the change of the cross sections with Q 2 . For comparison, we also 
show the photoproduction cross section in this figure (Q 2 = 0). From these plots, we can 
see that the shapes of the curves for the individual channel contributions do not change 
distinctly with Q 2 , and their relative importances also do not change. More important, 
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1 S'o channel always dominates over the other two channels, which is the same as that in 
the photoproduction processes. 

In Fig. 4, we show the ratio gl(i*p) / a T(l*p) as the function of z for Q 2 = lOGeV 2 and 
AOGeV 2 respectively. We can see that this ratio increases in all channels as Q 2 increases. 
In particular, for the S S[ 8 ^ channel, as z — > 1, the ratio o"l(7*p)/ct t (7*p) can reach a value 
over 4.4 for Q 2 = 40GeV 2 . This is a dictinct feature of the virtual photon with Q 2 » m 2 . 
But the total cross section <J to t{l*p) = 1l{i*p) + 0"t(t*p) coming from the 3 S[ 8 ^ channel 
decreases. The total cross section is alwayse dominated by the channel for all value of 

Q 2 . 

In Fig. 5, we show the cross sections of e + p — > e + 7] c + X at HERA. We plot the 
differential cross sections da/dz as a function of z for two cases for the integration region of 
Q 2 : (a) AGeV 2 < Q 2 < 80GeV 2 ; (b) 20GeV 2 < Q 2 < 80GeV 2 . 

From the above numerical results, we see that the contribution from the Sq subprocess 
dominates the r\ c electroproduction. The associated color-octet matrix element (Og c ( 1 S )) 
can be related to the color-octet matrix element {Og^ ( 3 Si)) by the heavy quark spin symme- 
try. The value of the later matrix element originates from the fitting to the large transverse 
momentum J/i/j production at the Fermilab Tevatron [|1^], which can be viewed as a reliable 
estimate. This is because the color-octet matrix element (0^( 3 Si)) dominantly contributes 
in the large px production [|,|]J§, where the initial and final state gluon radiation effects 
are much smaller 0. So, the r\ c electroproduction together with its photoproduction can 
provide an important test for the color-octet production mechanism. 

We note that the next to leading order QCD perturbative corrections are large for J/ip 
photoproduction at HERA energy ranges P JT7|JT8|] . These corrections change the normaliza- 
tion of the production rate by a factor of about 2. r\ c photoproduction and electroproduction 
may also be affected by these corrections. However, the main features in testing the color- 
singlet and color-octet mechanisms mentioned above for r\ c production will not be changed. 

For the experimental observation of 7] c , we suggest using all main decay channels of r\ c 
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such as rjTnr, i]'tttt, KKtt, tttttttt, KKtttt, KKKK, as well as 77 to detect r] c . The increase 
of luminosity at HERA offers the possibility of experimental investigation of eletrophotopro- 
duction of r\ c . In this connection, we note that photoproduction of the D s meson has been 
observed through the D s — ► 07r channel at HERA (T^|. According to our calculation, the 
photoproduction cross section of r] c is of the same order as the observed D s . The branching 
ratio of r] c — > KKtt or r] c — > 00 is comparable to that of Z) s — > 07r. Therefore, the experi- 
mental observation of photoproduction and leptoproduction of r\ c might be possible through 
i] c — > KKtt, 00 decay modes in the near future. 



IV. SUMMARY 

We have calculated in this paper the inelastic r\ c production in the electroproduction 
processes under the NRQCD factorization formalism. Same as the photoproduction of r\ c 



T2|| , the color-singlet contributions to the electroproduction also vanish in the leading order 



and the next to leading order in perturbative QCD. The dominant contribution comes from 
the color-octet 1 Sq 8 ^ subprocess, for which the associated color-octet matrix element is related 
to the direct production of J/ip at large transverse momentum at the Fermilab Tevatron. 
The measurement of this process at HERA can be viewed as another independent test for 
the color-octet production mechanism, and complement to the study of its photoproduction. 
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APPENDIX: 



In this appendix, we list the cross sections for the different partonic processes. 
For 1 S 3 Si and 1 Pi, the cross sections: 

F( 3 S{ 8 ^ = , m-?elaal 14 12 12£ 

1 1J 3(s + t + Qiy(s + u + 2Q 2 )i(t + u + Qiy(Q* + sy [ ^ v v 

+7Q 12 m + 22g 10 s 2 + 40g 10 st + 30g 10 sw + 16g 10 P + 18Q 10 iu + 5Q w u 2 
+ 18Q 8 s 3 + 55g 8 s 2 t + 50Q 8 s 2 u + 5ig 8 s£ 2 + 6AQ 8 stu + 2lQ 8 su 2 + 12Q 8 P 
+ 19Q 8 Pu + 8Q 8 tu 2 + Q 8 u 3 + 7Q 6 s 4 + 35Q 6 sH + 40Q 6 s 3 u + QlQ 6 s 2 P 
+8hQ & s 2 iu + 34Q 6 s 2 u 2 + 3ig 6 sP + 5QQ 6 sPu + 27Q 6 siu 2 + AQ 6 su 3 + QQH 4 

+ng 6 Fw + 7g 6 £ 2 w 2 + Q 6 iu 3 + g 4 s 5 + ioqHH + i5Q 4 s% + 34g 4 s 3 t 2 

+5ig 4 s 3 M + 26g 4 s 3 w 2 + 27Q 4 s 2 t 3 + 58Q 4 s 2 Pu + 33Q 4 sHu 2 + 6Q 4 s 2 u 3 

+8g 4 st 4 + 2ig 4 st 3 w + i7g 4 st 2 w 2 + 3g 4 sfe 3 + 2g 4 t 5 + 4g 4 t 4 « + 3g 4 £ 3 « 2 

+Q 4 i 2 u 3 + Q 2 sH + 2g 2 s 5 w + 9Q 2 s 4 P + 13g 2 s 4 M + 9g 2 sV + 9g 2 s 3 t 3 

+24g 2 s 3 t 2 « + i7g 2 s 3 fe 2 + 4g 2 s 3 # + 2g 2 s 2 t 4 + iiqWu + i3g 2 *W 

+3g 2 s 2 M 3 + 2g 2 st 4 w + 4g 2 st 3 w 2 + 2Q 2 st 2 u 3 + sH 2 + sHu + s 5 u 2 + s 4 i 3 
+3sH 2 u + 3s 4 M 2 + s 4 u 3 + s 3 i 3 u + 3s 3 Pu 2 + s 3 iu 3 + s 2 i 3 u 2 + s 2 t 2 u 3 }. (Al) 



PC^) = ^ ; : ; ^ ' TZ", ^ , l , ^tW 6 + 45g 44 s + 4 9 g 44 t 



i c (8)^ 24(47r) 3 e;W 2 

(g 2 + s + t) 2 (2g 2 + s + u) 2 {Q 2 + s) 2 (g 2 + t + u)H 

+27Q 14 u + 9QQ 12 s 2 + 17AQ 12 st + 123Q 12 su + 9AQ 12 t 2 + 136g 12 M + 33g 12 w 2 
+ 114g 10 s 3 + 256g 10 s 2 £+ 234g 10 s 2 w + 250g 10 sP + 430g 10 stw + 135g 10 sw 2 

+8og 10 £ 3 + 202g 10 Pw + i45g 10 M 2 + 2ig 10 « 3 + 82g 8 s 4 + 206g 8 s 3 £ 

+242g 8 s 3 w + 273Q 8 s 2 t 2 + 538g 8 s 2 tw + 222g 8 s 2 w 2 + 172g 8 st 3 + 462g 8 st 2 w 

+406g 8 sM 2 + 77g 8 sw 3 + 3ig 8 £ 4 + ii2g 8 t 3 « + i45g 8 t 2 « 2 + 7og 8 t« 3 + 7g 8 « 4 

+36g 6 s 5 + 98g 6 s 4 t + 148g 6 s 4 w + 157g 6 s 3 £ 2 + 352g 6 s 3 M + 188g 6 s 3 w 2 
+ 153g 6 s 2 t 3 + 413Q 6 s 2 Pu + 420g 6 s 2 M 2 + 108g 6 s 2 w 3 + 63g 6 s£ 4 + 194g 6 st 3 w 
+275g 6 s£ 2 w 2 + 172g 6 sfe 3 + 23g 6 sw 4 + 5Q 6 i 5 + 2ig 6 t 4 w + 39g 6 t 3 w 2 
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+35Q 6 t 2 w 3 + 13Q 6 tw 4 + Q% 5 + 9Q 4 s 6 + 2QQ 4 s 5 t + 54Q 4 s 5 u + 45Q 4 s 4 P 
+132Q 4 s 4 iu + 87Q 4 s 4 u 2 + 60Q 4 s 3 t 3 + 191Q 4 sH 2 u + 208Q 4 s 3 iu 2 + 72Q 4 s 3 u 3 
+43Q 4 s 2 t 4 + 134Q 4 s 2 t 3 w + 18ig 4 sW + 140Q 4 s 2 tw 3 + 27Q 4 s 2 u 4 + 10Q 4 st 5 
+33Q 4 st 4 w + 50Q 4 si 3 u 2 + 53Q 4 si 2 u 3 + 28Q 4 siu 4 + 3Q 4 su 5 + Q 2 s 7 + 3Q 2 sH 
+ llQ 2 s 6 u + 5QW + 2QQ 2 sHu + 21Q 2 s 5 u 2 + 9Q 2 sH 3 + AlQ 2 sH 2 u 
+hhQ 2 sHu 2 + 23Q 2 s 4 m 3 + lig 2 s 3 t 4 + 34Q 2 s 3 Pm + 57g 2 s 3 t 2 w 2 + 44g 2 s 3 tw 3 
+13QW + 5Q 2 sH 5 + 13Q 2 s 2 t 4 w + 23Q 2 s 2 Fm 2 + 21Q 2 s 2 t 2 w 3 + l7Q 2 s 2 tu A 
+3Q 2 s 2 u 5 + fu + 2sHu + 2s 6 m 2 + 3s 5 Pu + 6§Hu 2 + 3s 5 u 3 + 2sH 3 u + 6s 4 Pu 2 
+6s 4 iu 3 + 2s 4 u 4 + s 3 i 4 u + 2s 3 t 3 w 2 + 3s 3 Pw 3 + 2s 3 tu 4 + s 3 u 5 ]. (A2) 

320(4ir) 3 e 2 c aa 2 

3(2m c ) 2 (g 2 + I + t)3(2g 2 + I + w) 4 (g 2 + t + u)\Q* + s) 2 

[48Q 22 + 336g 20 s + 336g 20 t + 240g 20 w + 1048g 18 s 2 + 1900g 18 st + 1552g 18 sw 

+922g 18 £ 2 + 1548g 18 M + 520g 18 w 2 + 1920g 16 s 3 + 471QQ w sH + U32Q w s 2 u 

+4304g 16 st 2 + 7960g 16 stw + 3088g 16 s« 2 + i4iog 16 P + 3700g 16 £ 2 w 

+3044g 16 tw 2 + 640g 16 w 3 + 229ig 14 s 4 + 6795g 14 s 3 t + 7356g 14 s 3 w 

+8702g 14 s 2 t 2 + 17713g 14 s 2 M + 7994g 14 s 2 w 2 + 5498g 14 sP + 15470g 14 st 2 w 

+14125g 14 sM 2 + 3468g 14 sw 3 + 1518g 14 t 4 + 4898g 14 Fw + 6220g 14 £V 

+3319g 14 fe 3 + 49ig 14 w 4 + 1863g 12 s 5 + 6339g 12 s 4 £ + 7853g 12 s 4 w 

+10120g 1 W + 22468g 12 s 3 M + 11838g 12 s 3 w 2 + 9156g 12 s 2 t 3 + 27460g 12 s 2 Pw 

+27824g 12 s 2 £w 2 + 8050g 12 s 2 w 3 + 4972g 12 s£ 4 + 16826g 12 sF« + 22970g 12 sPw 2 

+ 13768g 12 sM 3 + 241ig 12 sw 4 + 1334g 12 t 5 + 4676g 12 £ 4 « + 7100g 12 Pw 2 

+5690g 12 t 2 w 3 + 2177g 12 M 4 + 24ig 12 w 5 + 1047g 10 s 6 + 4042g 10 s 5 t 

+5628g 10 s 5 w + 7609g 10 s 4 t 2 + 18020g 10 s 4 tw + 11052g 10 s 4 w 2 

+8639g 10 s 3 £ 3 + 27282g 10 s 3 £ 2 w + 30485g 10 s 3 M 2 + 10452g 10 s 3 w 3 + Q717Q 10 s 2 t 4 

+23967g 10 s 2 £ 3 w + 35078g 10 s 2 Pw 2 + 23615g 10 s 2 M 3 + 4955g 10 s 2 « 4 
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+3568Q 10 sF + 13176Q 10 st 4 M + 20989Q 10 sPu 2 + 18212Q l0 sPu 3 + 7977Q w stu 4 
+ 1064Q 10 sm 5 + 902Q 10 t 6 + 3520g 10 £ 5 w + 5917Q 1( W + 5581Q 10 ?m 3 
+3079Q 10 ?m 4 + 869Q 10 fe 5 + 74Q 10 u 6 + 402Q 8 s 7 + 1807Q 8 s 6 t 
+2736Q 8 s 6 m + 3946Q 8 s 5 t 2 + 95Q8Q 8 s 5 iu + Q747Q 8 s 5 u 2 + 5199QW 
+16875QWu + 20452Q 8 s 4 M 2 + 8305Q 8 s 4 m 3 + 4895g 8 s 3 t 4 + 18541QWu 
+28923Q 8 s 3 Pm 2 + 21762g 8 s 3 M 3 + 5520g 8 s 3 w 4 + 3687Q 8 s 2 t 5 + 14625Q 8 s 2 t 4 w 
+24778g 8 s 2 £ 3 w 2 + 23241Q 8 s 2 t 2 w 3 + 11647g 8 s 2 fe 4 + 1906Q 8 s 2 m 5 + 1848Q 8 st 6 
+7722g 8 st 5 w + 13765Q 8 st 4 w 2 + 13745g 8 sPw 3 + 8329g 8 st 2 w 4 + 2770Q 8 sM 5 
+291Q 8 su 6 + 398g 8 f + 1804Q 8 t 6 w + 3523g 8 t 5 w 2 + 3855Q 8 t 4 w 3 + 2547g 8 t 3 w 4 
+998g 8 t 2 w 5 + 202Q 8 M 6 + 13g 8 w 7 + 10ig 6 s 8 + 5Q7Q 6 s 7 i + 887Q 6 s 7 u 
+ 1458g 6 s 6 P + 3403g 6 s 6 tw + 2693g 6 s 6 w 2 + 2141QW + 6828QWu 
+8693g 6 s 5 tw 2 + 4139g 6 s 5 w 3 + 2154g 6 s 4 t 4 + 8625g 6 s 4 Pw + 14103g Ww 2 
+ 11713g 6 s 4 M 3 + 3600g 6 s 4 w 4 + 1868g 6 s 3 £ 5 + 8205g 6 s 3 t 4 w + 14993g 6 s 3 £ 3 w 2 
+ 15220g 6 s 3 £ 2 w 3 + 8678g 6 s 3 tw 4 + 1769g 6 s 3 w 5 + 1337g 6 s 2 t 6 + 6118g 6 s 2 t 5 w 
+ 11895g 6 s 2 t 4 w 2 + 1281ig 6 s 2 t 3 w 3 + 8502g 6 s 2 t 2 w 4 + 3332g 6 s 2 M 5 + 445g 6 s 2 w 6 
+570g 6 st 7 + 2792g 6 st 6 w + 5910g 6 st 5 w 2 + 6989g 6 st 4 w 3 + 5004g 6 st 3 w 4 
+2210g 6 st 2 w 5 + 550g 6 sM 6 + 45g 6 sw 7 + 98g 6 t 8 + hUQH 7 u + 1187g 6 t 6 w 2 
+ 1572g 6 t 5 w 3 + 1287g 6 tV + 6A6QH 3 u 5 + 183Q 6 t 2 u 6 + 24g 6 fe 7 + Q 6 u 8 
+ 15g 4 s 9 + 12ig 4 s 8 £ + 182g 4 s 8 w + 379QW + 79ig 4 s 7 fe + 677g 4 s 7 w 2 
+620g 4 s 6 P + 1823g 4 s 6 £ 2 w + 2317g 4 s 6 M 2 + 1265g 4 s 6 w 3 + 623g 4 s 5 t 4 
+2536g 4 s 5 t 3 w + 4166g 4 s 5 tV + 3725g 4 s 5 tw 3 + 1376g 4 s 5 w 4 
+498g 4 s 4 t 5 + 2545Q 4 s 4 t 4 u + 4986g 4 s 4 t 3 « 2 + 5449gWw 3 
+3502g 4 s 4 tw 4 + 896g 4 s 4 w 5 + 393g 4 s 3 £ 6 + 2150g 4 s 3 t 5 w + 4713g 4 s 3 t 4 w 2 
+5603g 4 s 3 t 3 w 3 + 4117g 4 s 3 t 2 w 4 + 188ig 4 s 3 tw 5 + 329g 4 s 3 w 6 + 237g 4 s 2 t 7 
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+1343Q 4 s¥w + 3229Q 4 s 2 t 5 w 2 + 4283Q 4 s 2 t 4 w 3 + 3395Q 4 s 2 Pw 4 
+1683Q Ww 5 + 519Q 4 s 2 tw 6 + 57Q 4 s 2 u 7 + 7QQ 4 st 8 + 460Q 4 st 7 w 
+ 1205Q 4 st 6 w 2 + 1804Q 4 st 5 w 3 + 1670Q 4 st 4 w 4 + 955Q 4 st 3 w 5 
+320Q 4 si 2 u 6 + 55Q 4 siu 7 + 3Q 4 su 8 + 10Q 4 t 9 + Q0QH 8 u + 1Q1QH 7 u 2 
+255Q 4 t 6 u 3 + 261QH 5 u 4 + 174Q 4 t 4 u 5 + 71Q 4 Pw 6 + 15Q 4 Pu 7 + QHu 8 
+Q 2 s 1 + 16Q 2 s 9 i + 21Q 2 s 9 u + 63Q W + llOQ 2 s 8 iu + 97Q 2 s 8 u 2 
+ 118QW + 302QWu + 359Q 2 s 7 iu 2 + 217Q 2 s 7 u 3 + 125QW 
+466Q 2 s 6 t 3 w + 715Q 2 s 6 tV + Q59Q 2 s 6 iu 3 + 286Q 2 s 6 u 4 + 8AQ 2 sH 5 
+m\Q 2 sH% + 910QWw 2 + 1032Q 2 s 5 Pw 3 + 735Q 2 s 5 tw 4 + 235Q 2 s 5 w 5 
+47Q 2 s 4 t 6 + 3UQ 2 sH 5 u + 876Q 2 sW + 1148Q 2 sW + 949Q 2 s 4 t 2 w 4 
+501Q 2 s 4 tw 5 + l\7Q 2 s% % + 26Q 2 s 3 f + 206Q 2 s 3 t 6 w + 640Q 2 s 3 t 5 w 2 
+999Q 2 s 3 £ 4 m 3 + 916QWw 4 + 528QWm 5 + 196Q 2 s 3 iu 6 + 31Q 2 s 3 u 7 
+8Q 2 s 2 t 8 + 78Q 2 s 2 t 7 u + 289Q 2 s 2 t 6 u 2 + 5A2Q 2 s 2 t 5 u 3 + 599Q 2 s W 
+402Q 2 s 2 i 3 u 5 + 161Q 2 s 2 i 2 u 6 + 38Q 2 s 2 iu 7 + 3Q 2 s 2 w 8 + 10Q 2 si 8 u + 52Q 2 si 7 u 2 
+120Q 2 si 6 u 3 + lQ2Q 2 si 5 u 4 + 138Q 2 st 4 w 5 + 72Q 2 st 3 w 6 + 20Q 2 sPu 7 + 2Q 2 siu 8 
+s 1 0t + s 4 0m + 5s 9 P + 7 sHu + 6s 9 u 2 + lls 8 t 3 + 2As 8 t 2 u + 2hs 8 iu 2 
+ 16s 8 m 3 + Uft 4 + A3s 7 i 3 u + 57s 7 i 2 u 2 + 5ls 7 iu 3 + 25s 7 w 4 + lls 6 t 5 + 46s 6 t 4 w 
+78s 6 Pw 2 + 84s 6 t 2 w 3 + 64s 6 M 4 + 25s 6 w 5 + 5s 5 t 6 + 30sH 5 u + 70s 5 t 4 w 2 
+92s 5 t 3 w 3 + 84s 5 tV + 51s 5 tu 5 + 16s 5 w 6 + sH 7 + llsH 6 u + 38s 4 t 5 w 2 + 70s 4 t 4 w 3 
+78s 4 Pw 4 + 57s 4 t 2 w 5 + 25sHu 6 + Qs 4 u 7 + 2s 3 i 7 u + lls 3 t 6 u 2 + 30s 3 Fw 3 
+46s 3 tV + 43s 3 Fw 5 + 24s 3 £ 2 m 6 + 7sHu 7 + s 3 u 8 + s 2 t 7 w 2 + 5s 2 i 6 u 3 + lls 2 i 5 u 4 
+UsH 4 u 5 + llsW + 5sW + s 2 tu 8 ]. (A3) 
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Figure Captions 



FIG. 1. The Feynman diagrams for the production of i] c at the leading order (a) and 
the next to leading order (b)and(c) 7*g fusion processes. All of these diagrams are for the 
color-octet processes. 

FIG. 2. The cross section da/dQ 2 for the process e + p^e + r] c + Xm the forward 
direction (z ~ 1) as a function of Q 2 . where v /s7^ = 300GeV. 

FIG. 3. The differential cross sections da/dz for the process 7*+p — * i] c +X at the HERA 
as a function of z = E^/E*, where W rp = WOGeV; (a) Q 2 = OGeV 2 ; (b) Q 2 = AGeV 2 ; 
(c) Q 2 = lOGeV 2 ; (d) Q 2 = 40GeV 2 . The curves denote contributions from the color-octet 
1 S , 1 Pi, and 3 S'i respectively. 

FIG. 4. The ratio a L (Yp) I '°t(i*p) as the function of z. (a) Q 2 = lOGeV 2 . (b) Q 2 = 
40GeV 2 . 

FIG. 5. The cross section da/dz for e + p — > e + r] c + X at HERA as a fuction of where 
fsT p = 300GeV. (a) 4GeV 2 < Q 2 < 80GeV 2 , (b) 20Ge\/ 2 < Q 2 < 80GeV 2 . 
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Y+P->ti c +X (W yp =100GeV) 




Fig.3 
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e+p->e+r| c +X (E mc =300GeV) 
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